Can a 3+2 Oscillation Model Explain the NuTeV Electroweak Results? 



J. S. MeQ 

Carnegie Mellon University 
Research Experience for Undergraduates at Columbia University 



in 
o 
o 

(N 



O ■ 

in ■ 
o ; 

■ 

Oh 

CD 



X 



J. M. Conrad, M. Sorel, G. P. Zeller 
Columbia University 
(Dated: February 7, 2008) 

The sin^ dw result from NuTeV falls three standard deviations from the value determined by global 
electroweak fits. It has been suggested that one possible explanation for this result could be the 
oscillation of electron neutrinos in the NuTeV beam to sterile neutrinos. This article examines several 
cases of masses and mixings for 3+2 neutrino oscillation models which fit the current oscillation data 
at 99% CL. We conclude that electron to sterile neutrino oscillations can account for only up to a 
third of a standard deviation between the NuTeV determination of sin^ 9w and the standard model. 
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NuTeV precisely determines the electroweak mixing 
angle through the measurement of deep inelastic muon 
neutrino and antineutrino interactions. Although the fi- 
nal value of sin^ 9\y obtained by this experiment agrees 
with previous neutrino-based measurements, the result 
is anomalously high when compared to the value from 
global electroweak fits to other data. The NuTeV result, 
0.2277 ± 0.0013(stat) ± 0.0009(sys) p], is roughly three 
standard deviations above the standard model value of 
0.2227 ±0.0004 ,2]. 

Giunti et al. jj] have considered neutrino oscillations as 
a possible explanation for the NuTeV results, suggesting 
that if electron neutrinos in the NuTeV beam were os- 
cillating into sterile neutrinos, this could effectively lead 
to the NuTeV observation. Their paper demonstrated 
that a 34-1 (three active and one sterile) neutrino model 
would require very large mixings to the sterile neutrino. 
Such large mixings are now known to be inconsistent with 
present oscillation limits. In addition, the proposed oscil- 
lations are too large to be consistent with the direct mea- 
surement of the electron neutrino content in the NuTeV 
beam In this paper, we extend the idea in Q to oscil- 
lation models with two sterile neutrinos, i.e. 3-1-2 (three 
active and two sterile) neutrino models. For a reveiw of 
these models and their motivation, see Reference 



I. THE NUTEV DETECTOR AND ANALYSIS 

The design of the NuTeV experiment is described in 
detail in Reference 0. This experiment used a high 
energy 800 GeV proton beam, taking data in neutrino 
and antineutrino modes separately. The NuTeV detector 
was located 1450 m downstream from the proton target, 
and consisted of a steel-scintillator target followed by a 
toroid spectrometer. Two types of interactions can occur: 



charged current events (CC), which proceed by ex- 
change, and neutral current (NC), which proceed by 
exchange. Both interactions produce a hadron shower of 
particles in the calorimeter. For NC events, the shower is 
accompanied by an undetectable final state neutrino. For 
CC events, there is instead a muon which can be tracked 
through the calorimeter and the toroid spectrometer. To 
lowest order in both QCD and electroweak theory, the 
ratio of NC to CC rates in neutrino and antineutrino 
scattering relates directly to sin^ 9w : 
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where R'^''^ is the ratio of NC to CC total cross sections 
and r is the ratio of muon neutrino to antineutrino CC 
total cross sections. 

To extract a value of sin^ 6w from the data, NuTeV 
does not measure total cross section ratios, but rather 
measures experimental ratios of NC to CC candidate 
events. NuTeV differentiates NC and CC interactions 
simply by the measured length of the event NC 
interactions strictly produce hadronic showers and ap- 
pear as short events in the detector. Longer events are 
likely to be extended by virtue of containing a muon, 
and thus are identified as CC events. The total number 
of short and long events (iV^p and N^^p, respectively) 
are measured and from them, the experimental ratio, 
Rexp = Nf^p/N^^p is determined in both the neutrino 
and antineutrino data. These ratios include the effects of 
experimental cuts, cross-talk between candidates in the 
numerator and denominator, final state effects, and non- 
muon neutrino backgrounds. 

The second largest background to iV^p, accounting for 
~ 5% of short events in neutrino mode and 6% in an- 
tineutrino mode, results from electron neutrino contam- 
ination in the beam, the dominant source of which are 
Kf^ decays. The electron neutrino background is deter- 
mined using beam Monte Carlo tuned to the neutrinos 
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observed from Kj^2 decays. This beam prediction is then 
checked against a direct measurement of the electron neu- 
trino content in the NuTeV data 5]. The predicted and 
measured electron neutrinos are found to agree: the ra- 
tio of measured to Monte Carlo predicted Vf. events is 
1.05 ± 0.03 in the neutrino beam and 1.01 ± 0.04 in the 
antineutrino beam 

Because of their event topology, electron neutrino in- 
teractions all appear as short events in the NuTeV detec- 
tor. Their contribution (in addition to other corrections 
which are not shown here) is explicitly included as a mod- 
ification to the number of short events appearing in the 
numerator of the predicted experimental ratio: 
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in both neutrino and antineutrino modes. Any overes- 
timate of the electron neutrino contribution (for exam- 
ple, that would result from neglecting — > Vg oscilla- 
tions) would lead to an overestimate of the predicted ra- 
tio, Reip , and hence a larger measured value of NuTeV 
sin^ 9w- Importantly, any adjustment to the electron 
neutrino flux needed to reduce the NuTeV sin^ 9w value 
and bring the result into better agreement with expecta- 
tion must additionally satisfy the direct constraint from 
the NuTeV data itself ^J. Here, we consider several such 
possibilities. 



II. NEUTRINO OSCILLATIONS 

To gain an understanding of the potential impact of 
electron neutrino oscillations in NuTeV, first consider the 
approximation of a simple two-flavor (i^e — > J^s) oscilla- 
tion probability: 



P = sin^ 2esi-a^{l.27Am'^L/E). 



(4) 



The fundamental parameters describing the oscillation 
are sin^ 26*, the mixing between the flavors, and Am^, 
the squared mass difference between the neutrinos. The 
experimental parameters are L, the baseline, and E, the 
incident neutrino energy. In general, oscillations become 
observable when Arn^L/E '--^ 1 or larger. The high beam 
energy and short baseline of NuTeV lead to a small value 
of L/E, therefore requiring a large value of Am^ (of a 
few eV^ or greater) to compensate. 

In the few eV^ range of Am^, the Bugey reactor ex- 
periment sets the best limit on the mixing angle for 
disappearance 0. Reference shows that, in a 3-1-1 
model, the NuTeV result implies Am^ ~ 10 eV^ and 
~ sin^ 20 = 0.4, which would have produced a clear sig- 
nal in Bugey; therefore, this solution is directly excluded. 
The situation for 3-1-2 models, however, could differ. In 
such a model, the probability for electron neutrino disap- 
pearance at high Am^, with mixing matrix parameters 
denoted as U rather than mixing angles, is given by: 
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FIG. 1: 90% (dark grey) and 99% (light grey) CL allowed 
regions in (Am|i, (7e5)-space for CP-conserving (3-1-2) models. 
The stars labeled A, B, C, and D indicate the four models 
evaluated for their impact on the NuTeV sin^ 9w analysis. 



where Xij = 1.27 • Amfj{L/E). Such 3-1-2 models can 
fit the world's oscillation data with many different val- 
ues of mass splittings and mixing parameters fl. Those 
which provide a good description of the data tend to have 
Amj^ - 1 eV^ and Am^i > 10 eV^. 

For NuTeV, which has a small L/E, in the case where 
Amli << 10 eV^, (which is the 3-1-2 best fit case), and 
taking 2:41 and X54 — X51, the oscillation probability 
simplifies to: 



P-4[/f5(l-C/f5)sin2x5i 
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analogous to Equation 0] Figure Q shows the 90% and 
99% CL allowed regions in the Amli and Ue5 mixing 
space, as determined using the methods described in 
Reference @. The short-baseline data yielding the pa- 
rameter space relevant for this study include Bugey 
CCFR84n3, CDHS [13 , CHOOZ O, KARMEN2 O, 
LSND TilTand NOMAD 



III. IMPACT OF 3+2 MODELS ON THE 
NUTEV ELECTROWEAK RESULTS 

Given the increase in possible parameter space inher- 
ent in 3+2 models, we ask whether there is an oscillated 

flux which can account for the NuTeV electroweak re- 
sults. To answer this question and demonstrate the im- 
pact of a 3+2 model neutrino oscillations on the NuTeV 
data, four representative points are selected within the 
allowed region as indicated in Figure Q These points 
are the best fit model (A), a high mixing model with 
lower mass (B), a high mixing model with higher mass 
(C), and a high mass model (D). The best fit point (A) 
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TABLE I: Comparison of NuTeV electroweak results assum- 
ing no Ve oscillations (default) to the results assuming three 
CP-conserving 3+2 oscillation models (B,C,D in Fig. 0. In 
all cases, the same event selection criteria as in Q is applied. 
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FIG. 2: Unoscillated NuTeV i/e flux prediction in each mode. 
Inlay shows the ratio of the predicted oscillated/unoscillated 
Ve fluxes as a function of energy for Model D. The change is 
within the errors of the NuTeV measurement Q. 

was eventually omitted from the study because it rep- 
resented such a small correction to the unoscillated flux 



(as expected from Equation 0) that it had a negligible 
impact on NuTeV sin^ 6w . 

For each set of possible parameters, a j/g survival prob- 
ability is calculated as a function of energy. This prob- 
ability is then used to correct the estimated NuTeV i^e 
and Ve fluxes. The resulting difference is very small in 
all cases, being largest for the high mass Model D (Fig- 
ure [5J. The total integrated i/^ flux prediction changes 
by 0.2% (Model B), 0.8% (Model C), and 1.8% (Model 
D), hence satisfying the NuTeV z^e data constraint |^. 
Based on each of the "oscillated" i^e and fluxes, the 
Reip predictions are then recalculated in both neutrino 
and antineutrino modes, and a new value of sin^ 9yy is 
extracted. In Table HJ we report the deviation of the 
measured R^^p and sin^ 9\y values for the three mod- 
els: B, C, and D. The magnitude and sign of the shift 
indicates the expectation if the NuTeV data had been an- 
alyzed using the oscillated electron neutrino fluxes. All 
of the adjusted fluxes move the NuTeV results into bet- 
ter agreement with the standard model, by construction. 
Model D provides the largest impact: shifting R^^p and 
into better agreement with expectation by 0.4 a and 

0.3 a, respectively. Model D reduces the NuTeV sin^ 9w 
discrepancy with the standard model from S.Oct to 2.7 a. 

We have studied three points representative of extreme 
masses and mixings within the allowed region of 3+2 
models, and conclude that Ve Vs oscillations in this 
model do not yield a significant impact on NuTeV's elec- 
troweak results. The largest shift in sin^ 9w is created 
by a high mass model (e.g. model D), but even such 
a high mass model would only affect the NuTeV value 
of sin^ 9w by roughly 0.3 ct. Therefore, a 3+2 model 
with Vf, — > Vg oscillations cannot explain the NuTeV elec- 
troweak results by itself. 
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